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The (-( Correlator Is Time Dependent 
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We comment on the recent arguments by Senatore and Zaldarriaga that loop corrections to the 
(,-(, correlator cannot grow with time after first horizon crossing. We first emphasize the need to 
search for such secular dependence in corrections whose in-out matrix elements are infrared singular 
on an infinite spatial manifold. Then we give examples of such time dependence from pure quantum 
gravity and from scalar potential models. Finally, we point out that this time dependence arises 
from inflationary particle production and is therefore unlikely to endanger the preservation of super- 
horizon correlations as a record of inflation. 

PACS numbers: 98.80.Cq, 04.62.+V 



INTRODUCTION 



Perhaps the most commonly quoted result for models 
of inflation is the curvature power spectrum [l[ , 
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(The field 1Z is defined by stripping the derivatives from 
the 3-curvature in the co-moving frame for which the mo- 
mentum flux vanishes [lj.) These predictions are made in 
the context of perturbation theory about a homogeneous, 
isotropic and spatially flat geometry 

o, H 

ds 2 = -dt 2 + a 2 (t)dx ■ dx =► H = - , e = (2) 

a H z 

The time of first horizon crossing is tk such that k = 
H(tk)a(tk), after which A^_(k,t) becomes nearly con- 
stant and one drops the argument t. The tree order result 
for typical single-scalar inflation models is 
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Theorists are eager to predict A^(fc) because its value 
for cosmological wave lengths can be reconstructed from 
observations of anisotropies in the cosmic microwave ra- 
diation and from large scale structure surveys [3] , 
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-0.037±0.012 



(4) 

This connection between quantum gravity and cosmolog- 
ical observation represents one of the great triumphs of 
inflation theory [ij , and accords expression (Q| the status 
of the first quantum gravitational data ever obtained. 

Tree order results such as ^ derive from the lin- 
earized mode functions. There can be important contri- 
butions from times before tk [5|> but the mode functions 
become constant afterwards because the restoring force 
k 2 /a 2 (t) redshifts away while the friction term remains 



large. Quantum loop effects can induce late time de- 
pendence by coupling mode k to changes in the vacuum 
energy from the quantum fluctuations of other modes. 
A theorem by Weinberg limits this time dependence to 
powers of the "infrared logarithm", ln[a(i)/a(ifc)] 
No one disputes this bound, the issue is its saturation. 

In his first paper on the subject Weinberg considered 
two one loop processes which seemed to contribute in- 
frared logarithms [a]: 

• Section V gave a qualitative treatment of self- 
interactions within the gravity-inflaton system, cul- 
minating in equation (41); and 

• Section VII gave a computation of the contribution 
from Af free, massless, minimally coupled scalars, 
culminating in equation (71). 

Although other work has produced similar results 0, [l(| , 
Senatore and Zaldarriaga have argued that there can- 
not be any infrared logarithms from Weinberg's second 
source 



ll| . We agree 



indeed, this follows from a sim- 
However, 
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pie rule for counting infrared logarithms 
we do not accept the subsequent conclusion by Senatore 
and Zaldarriaga that A 2 z (k,t) is free of infrared loga- 
rithms from any source and to all orders. (A recent pa- 
per by Giddings and Sloth also disputes their conclu- 
sion [13].) The purpose of this paper is to show that in- 
frared logarithms arise at one loop from self-interactions 
of the gravity-inflaton system (as in Weinberg's first ex- 
ample) and at two loops from massless, minimally cou- 
pled scalars with a quartic potential. 

In section 2 we summarize the Lagrangian. Section 3 
describes a simple rule for counting the maximum num- 
ber of infrared logarithms which can derive from a given 
interaction [12J. In section 4 we compute a one loop ef- 
fect from self-interactions of the gravity-inflaton system. 
Section 5 gives a two loop effect from the potential of 
a massless, minimally coupled scalar. Our conclusions 
comprise the final section. 
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GAUGE-FIXED LAGRANGIAN 

The model we consider consists of three fields: the 
spacelike, D-dimcnsional metric g M!/ ; the scalar inflaton 
tp whose slow roll down its potential V((p) drives inflation; 
and a spectator scalar a which is centered at the <7q = 
minimum (of zero) of its massless potential U(a). The 
Lagrangian is, 
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U(a) 



(5) 



where R is the .D-dimensional Ricci scalar and a comma 
denotes ordinary differentiation. 

We decompose g M „ into lapse, shift and spatial metric 
according to Arnowitt, Deser and Misner (ADM) [lij ]. 

g^dx^dx" = -N 2 dt 2 +g ij {dx i -N i dt){dx j -N j dt) . (6) 

ADM long ago showed that the Lagrangian has a very 
simple dependence upon the lapse |l4j . 



C = ^Surface Tcrms^ — 
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The quantity A is a potential energy, 
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B 
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V(<p) + U(a) + ^9 ij ^, iiPt: 



(8) 

where R is the (D — l)-dimensional Ricci scalar formed 
from gij . The quantity B in ([7]) is a sort of kinetic energy, 



(9) 



B = {E\y-E %1 E ij -%KGyUp-<p, i N l ) +[a-ajN' 
where Eij/2N is the extrinsic curvature. 
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and a semi-colon denotes covariant differentiation. Vary- 
ing ([7]) with respect to TV produces an algebraic equation, 



(11) 



This gives the constrained Lagrangian a "virial" form 

V9 
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= ( Surface Terms 



8ttG 
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Further progress requires the use of perturbation the- 
ory. The nonzero background fields are g^ — a 2 (t)8ij 
and ip = ifo(t). The two nontrivial Einstein equations 
can be used to eliminate the background scalar, 



p-2) • (D-2) 
= — uz^r H » v ( ( Po) - 
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(13) 



Note that the background values of the potential and 
kinetic terms are equal, A = B = (D-2)[H+(D-1)H 2 }. 
Hence the background value of the lapse is unity. 
We fix time as Maldacena [l5| and Weinberg Q , 



G (t,x) = (p(t,x) - ip (t) = 



(14) 



The other (D — 1) conditions have to do with how we 
define the unimodular part of the metric gu, 
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We require g^ 



(15) 



hij to be transverse, 



Gi(t,x) = djgij(t,x) = djhij(t,x) = . (16) 

(Maldacena and Weinberg imposed transversality on the 
logarithm of gij ) The resulting Faddeev-Popov determi- 
nant depends only on hij, and is singular for e = 0. 

Of course no gauge can eliminate physical inflatons; 
with condition (fT4|) that degree of freedom resides in 
((t, x). Linearized gravitons are carried by hij(t,x), and 
spectator scalars are in <r(t,x). By contrast, the shift 
field N l (t,x) is a constrained variable which mediates 
interactions between the other fields. 

To reach a perturbative form we first employ (p~5)) to 
exhibit how the potential ((8|) depends on (, hij and cr, 



A = Aq-R+WttG 



[/(a)- 



-2C 



A (l + a) . 
(17) 



Here the spatial Ricci scalar is, 
e -2C r . 



R 



i?-2(L>-2)V 2 C-0D-2)(£>-3)C' fc C 



(18) 



where R = 0(h 2 ) is the Ricci scalar formed from g^ and 
V 2 = dig^dj is the covariant scalar Laplacian. At this 
stage we can also recognize that 1Z is just C, in D = 4 
dimensions and to linearized order [l[ , 



at,x)+o(exh,h< 



^, ~s a 2 (t) n (D-2 

The kinetic energy ^ can be expressed as, 
B = A + 2(D-2)H[(D-l)((-t k N k )-N k k 
+(D-2)(C-^ k N k )\(D-l)(C-^ k N k )-2N' 



(19) 



./,■ 



+(A fe fe ) 2 - E M E M - 8irG(&-<r :k N k ) 2 
Here we define N = N\ A, = g %0 N^ and E l} = ±[N i;j 



(20) 
(21) 
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The next step is to expand the volume part 



of the constrained Lagrangian in powers of a and (3, 
V5 a^-V 13 - 1 ^ 
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As Weinberg noted, the terms involving no derivatives of 
the gravity fields sum up to a total derivative [6j . Another 
important fact is that quadratic mixing between N l and 
£ can be eliminated with the covariant field redefinition, 

S k = N k + 9 ke d e ^ [^Jv 2 C - e(C-O^)] ■ (24) 

After much work the quadratic Lagrangians emerge, 

D-i , . _. , n _9±„ _ „ x 

d e S k 8 k S e j, (25) 

(D-2)ea D - > 
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^ = 64ttG ( ^ ^ ~~ a 2 9k hij dk hij ) ' 
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{ct 2 - -^dfcudfccr}. 



(26) 
(27) 
(28) 



Expression (1281) reveals cr to be a massless, minimally 
coupled scalar with unit normalization. Let us call its 
propagator iA(x;x'). From (|27|) . and relations f| 15til6|) . 
we see that the graviton propagator is proportional, 



(x; x') = 32itG 



D-2 



iA(x;x') 



(29) 

where = Sij — didj/S7 2 is the transverse projection 
operator. These relations are exact. Because e(t) = 
—H/H 2 is nearly constant during inflation, expression 
(|26p implies a similar relation for the £ propagator, 

The massless, minimally coupled scalar has a well- 
known infrared problem [16j | which we regulate by work- 
ing on T ^ 1 with radius L and then making the integral 
approximation for the mode sum 17] . 



iA(x;x') 



iD-l 



(2tt) 



x {e(t-t')u(t, k)u*(t', k) + 9(t'-t)u*(t, k)n(t', k)\. (31) 
The mode function for constant e is, 
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Constant e implies Ha e is also constant and hence, 

H(t k ) 



D-4=0=e 



lim u(t, k) = C(e) x 



(33) 

where C(0) = 1 and we will use C(e) ~ 1 generally. 

Relations (IT91 , (j3"0")) and (|33l) allow a trivial derivation 
of the typical result ([3]) for the scalar power spectrum, 
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The tensor power spectrum is, 



A 2 (fc,i) 



2ir 2 



d 3 3 



n 



hij(t,x)hij(t, 0) 



n 



(35) 

Relations (|29p and (|33p offer a similarly straightforward 
derivation of the typical tree result for A^(fc, t) 0, 0, 



A 2 (fc,i) 



fc 3 1 fi 

= * x 32ttGx 2 x \u(t,k)\ 2 ^—GH 2 (t k ) . 

tree 2?H 7T 

(36) 

The 1/e enhancement of the scalar power spectrum with 
regard to the tensor one presumably explains why the 
scalar contribution has been detected but the tensor sig- 
nal has so far not been resolved. At 95% confidence the 
bound on their ratio at fco = 0.002 Mpc -1 



is 



A 2 , (fc ) 



< 0.22. 



(37) 



With the typical tree order results © and (l3l)l) . and the 
measured spectrum ^ , this bound implies an upper limit 
on the inflationary Hubble parameter, 



GH 2 (t ko )^-xrx A^fco) < 10- 10 . 
lb 



(38) 



We can also get a bound on e by combining the typical 
tree results © and |36j with §37}, 



(39) 



e(t ko ) « — < 0.014 



This is why it was justified to use C(e) ~ 1. We shall go 
further and approximate loop corrections by taking the de 
Sitter limit of e = once multiplicative factors of e have 
been removed from vertices and propagators. 

It remains to derive the relevant interactions by ex- 
panding the constrained Lagrangian (|23[) . These interac- 
tions are quite complicated, but most of them are pre- 
cluded by too many differentiated fields from contribut- 
ing the maximum number of infrared logarithms. If we 
want just the maximum possible number of infrared log- 
arithms then the number of interactions at any order 
becomes manageable. To study the lowest order effects 
of £ self-interactions it suffices to consider the minimal 
generalization of (|26|) . 



(D-2)e g _ x (p _ 1K J7 a 
c WttG 



{c 2 - ^d k Cd k (} . (40) 



To study the lowest order effects of the scalar potential 
we need only, 
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INFRARED LOGARITHMS 



Infrared logarithms are factors of ln[a(t)] which can 
contaminate loop corrections involving undifferentiated 



gravitons or massless, minimally coupled scalars. The 
oldest example is from 1982 [18j and consists of the 
coincidence limit of iA(x;x') on de Sitter background 
(e(t) = and H{t) = Hi). With full dimensional regu- 
larization, and L = (_ff/a/) _1 , the result is 17|. 
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(42) 
(43) 

(44) 



+ K^) + ^- i)+Vj(i)+0 (K)}- (45) 

This exhibits the fallacy of the argument Senatore and 
Zaldarriaga gave against infrared logarithms based on 
"making the integral dimensionless" [ll]]. That is the 
change of variables from fc to z — k/Hia(t) in passing 
from (|4"3"|) to fl"4")) . Had the lower limit been k = this 
would indeed have eliminated any time dependence, how- 
ever, the integral would have been infrared divergent. So 
we come to a crucial insight: infrared logarithms derive 
from diagrams that would be infrared divergent as in- out 
matrix elements on the spatial manifold R ^ 1 . Of course 
that is why they are called infrared logarithms. 

Note that any derivatives, with respect to space or 
time, would have eliminated the infrared logarithm in 
(|45|) . This observation has led to a very simple rule for 
inferring the maximum number of infrared logarithms 
which can come from a particular interaction // 
the interaction has a total of K undifferentiated gravitons 
and undifferentiated massless, minimally coupled scalars, 
after partial integration has been exhausted, then each 
correction involving two such interactions can produce as 
many as K infrared logarithms. 

This rule has been tested in a variety of explicit, fully 
dimensionally regulated and renormalized computations 
on de Sitter background using the Schwinger-Keldysh for- 
malism For a massless, minimally coupled scalar 
with a quartic self-interaction the rule correctly predicts 
the number of infrared logarithms in the expectation 
value of the stress tensor at one and two loop orders 



2C | , and in the one and two loop order self-mass-squared 
21 , ■ For scalar quantum electrodynamics the rule 
correctly predicts the infrared logarithms which are seen 
in the one loop vacuum polarization [23| and in the two 
loop scalar and electrodynamic field strengths, as well as 
the two loop expectation value of the stress tensor 24 1. 
For a Yukawa-coupled scalar the rule has been checked 
with the one loop fermion self-energy and with the 
expectation value of the coincident vertex at two loop 
order l 26]. The rule also gives the correct number of 





FIG. 1: One loop correction from cubic and quartic self- 
interactions of £. 



infrared logarithms in the one loop fermion self-energy 
from quantum gravity [2t| . 

Senatore and Zaldarriaga considered two models in de- 
tail. The first, given in their equation (11) and studied 
in section 3, consists of a Lagrangian containing only dif- 
ferentiated fields ll|. Any such interaction has K = 0, 
so the rule predicts no infrared logarithms, which is what 
they found. Their second model, given in their equation 
(75) and studied in section 4, was the same as Wein- 
berg's: gravity + inflaton + N massless, minimally cou- 
pled scalars with no potential ll|. This model shows 
infrared logarithms, both from its Q self-interactions and 
from interactions with undifferentiated ha fields. How- 



ever, Senatore and Zaldarriaga ignored those interactions 
because they give no parametric enhancement involving 
the potentially large number N. The scalar kinetic terms 
which were the object of their study consist of differ- 
entiated a fields with a complicated set of couplings to 
one £ field. Although there are certainly some K = 1 
terms present, cancellations make the resulting integrals 
infrared finite 0, [ll[ so the rule again predicts no infrared 
logarithms, and that is what they found. 

Let us now consider the two interactions (f4TJ|) and f4"Tj) 
given in the previous section. The general form of the £ 
self-interaction (|40| is ( K d(^dC, which is the same as for 
quantum gravity. We therefore expect that there should 
be a single infrared logarithm from a correction involv- 
ing two 3-point interactions, with K = 1, or from a single 
4-point interaction, with K — 2. These corrections cor- 
respond to the diagrams depicted in Fig. [I] and we will 
show in the next section that they indeed produce a single 
infrared logarithm. Supposing that the spectator poten- 
tial is quartic, we see that (|41[) contains an interaction of 
the form £ 2 c 4 . This has K = 6, so the rule predicts three 
infrared logarithms from a correction which involves one 
such interaction. The corresponding diagram is depicted 
in Fig. [2] and we will confirm that it does produce three 
infrared logarithms in the penultimate section. 



TIME DEP. FROM SELF-INTERACTIONS OF ( 

The two diagrams of Fig. Q] derive from expanding ex- 
pression (|40|) to cubic and quartic orders. The second 
of these diagrams is very similar to the computation fea- 
tured in section V of Weinberg's paper Q . As he noted, a 
field redefinition would make (|40l) free were it not for the 
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extra factor of e -2 ^ on the term with space derivatives. 
Things can be simplified by exploiting Weinberg's obser- 
vation that only the time derivative term contributes an 
infrared logarithm at one loop order [f|. We therefore 
make the field redefinition, 



D-l 



e~ C_i 



D-l 



■In 



D-l 
1 + — 



(46) 

and forget about the residual interactions involving spa- 
tial derivatives. 

The Z propagator is the same as the £ propagator (|3U1) . 
Hence the one £ loop correction to the C~C correlator is, 
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C(x)C(x') 
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C loop 



D-l\ 2 



) 2 (ft|^(x)ZV) 



'D-l\2 



+\z\x)Z 2 {x') + \z*{ X )Z{x')\q) , (47) 
jiA(x;x')iA(x';x') 



(D-2)eJ 

~[iA(x; x')] 2 + «A(x; x)iA(x; x')}- (48) 



+ 



We now set x M = (t, x) and x ,AI = (t, 0), and Fourier 
transform on x. It also makes sense to retain only the 
infrared logarithm terms because time independent con- 
tributions derive as well from derivative interactions of 
the same order as (|40|) which we have ignored. That is 
where the ultraviolet divergences reside, and they can be 
absorbed into BPHZ counterterms as usual. In the ab- 
sence of any condition for fixing the finite parts of those 
counterterms, the infrared logarithm terms are the only 
unambiguous prediction. The final result is, 



A 2 



7v 



(M) 



GH- 
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f 27GH 2 



0{G 2 H 4 )Y 
(49) 



We should mention that it is by no means clear what 
collection of fields represents the observed scalar power 
spectrum A 2 z (k 7 t). At tree order it suffices to use the 
correlator, and our result (l49l) is based on extending 
that correspondence to all orders. This definition affords 
a simple renormalization scheme because then the power 
spectrum is a noncoincident Green's function of a funda- 
mental field, and ordinary renormalization makes those 
finite. However, it is conceivable that the measured quan- 
tity is actually the correlator of some composite opera- 
tor such as (|19j) . in which case an additional, composite 
operator renormalization would be required. Nonlinear 
modifications of the observable can introduce additional 
infrared logarithms. For example, if the correct observ- 
able is the correlator of the field Z(t,x) defined in ex- 
pression (|46|) . then there are no infrared logarithms at 
one loop order. However, there does not seem any rea- 
son to suppose this, and even doing so would not prevent 
the appearance of infrared logarithms at higher orders. 



FIG. 2: Two loop correction from the interaction given in 
expression (|4Tj) . Solid lines represent £ and dashed lines rep- 
resent a. 



TIME DEP. FROM SPECTATOR POTENTIALS 

Let us assume U(a) = A(T 4 /4!. The diagram of Fig. [2] 
derives from the £ 2 cr 4 term of the interaction (|41[) . 



AC 



(D-l) 
48 



Xea D -\ 2 a 4 



(50) 



The Schwinger-Keldysh fl9| result for this diagram is, 



(Fig.© 
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l(D-2)e- 





-iA + _(x;y)iA + _(x';y)j 



d D y{iA ++ (x;y)iA ++ (x';y) 
iX(D-l) 



D-l 



iA(y;y)] 2 . (51) 



The propagator zA ++ (x;x') is the same mode sum as 
PTjl . whereas iA + _(x;x r ) has the same first line as (J^J) 
but the curly-bracketed expression on the second line is 
replaced by just u*(t, k)u(t' , k). 

We again take x^ = (t, x) and x ,fl = (t, 0), and Fourier 
transform on x, to obtain, 



d D-X. „-* 



i8n 2 (D-l)\G 2 



ds [a(s)] 



D-l 



(D-2) 2 e 

{[u(t,k)] 2 [u*{s,k)] 2 -[u*(t,k)] 2 [u( S ,k)] 2 yiA} 2 . (52) 



There is no point in retaining the divergent part of the 
coincident propagator (|45p . and continuing to work in 
D dimensions, unless we add the various counterterm 
diagrams. That exercise is identical to the published two 
loop computation of the expectation value of the a stress 
tensor [20;]. We will therefore retain only the leading 
infrared logarithm terms and take D = 4. 

Oscillations of the mode functions preclude a coherent 
effect before first horizon crossing. After horizon cross- 
ing one may take the long wavelength limit of the mode 
functions, 



i{t,k) 



Hence the curly-bracketed term of ([52^1 becomes, 
{[«(*, k)] 2 [u*( S , k)] 2 - [u*(t, k)] 2 [u( S , k)} 2 } 

k 2 



iH r 1 

""Fl 3 



1 



1 



a 3 (s) a 3 (t) 
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Putting everything together produces, 



-D-l 



leading log 

XG 2 H 5 f\ l a 3 (s) 



8n 2 ek 3 J t 
XG 2 H 4 



ds 



1 



a3(t) 



ln 2 [a( S )], (55) 



- .^..^{^["(t)] + subleadingj . (56) 
And multiplying by k 3 /2tt 2 gives the power spectrum, 



cr loops 7T6 



GH 2 r XGH 2 
I 48?r 3 



In' 



: (a) + 0(A 2 )}. (57) 



CONCLUSIONS 



We have shown that the £~C correlator acquires time 
dependent infrared log corrections, starting at one loop 
(|4"9")l from C self-interactions, and at two loops ((57)) from 
the quartic potential of a spectator scalar. There should 
also be infrared logarithms from dynamical gravitons, 
starting at two loops. 

The physical interpretation of the spectator effect (|57[) 
derives from the small increase in the vacuum energy as 
inflationary particle production pushes the a field up its 
potential U{a). The dimensionally regulated and fully 
renormalized result for this has been derived [23], but 
we can understand its effect on the £ mode functions by 
simply adding the Hartree approximation of (|50p to the 
free £ Lagrangian (|26|) in D — 4 dimensions, 



8ttG 



U 2 -\d k cd k Q 



ln 2 (a)C 2 } 



(58) 

After horizon crossing the associated mode equation is, 

3Hu: 



3XGH* 2 



XGH 2 , 3 
In 3 



(a)} . (59) 



Inserting the quantum corrected mode function in ex- 
pression ([Ml) gives precisely our result (|5T[) . It seems 
likely that a similar explanation can be given for the ef- 
fects from £ self-interactions, and from interactions with 
gravitons. 

Such effects must be present or else there is something 
seriously wrong with our understanding of how gravity 
responds to quantum fluctuations. That they would even 
be questioned is a tribute to how firmly cosmologists have 
come to believe in the time independence of £(t, k) after 
horizon crossing. Of course we appreciate the wonder 
of preserving a memory of conditions from inflation, but 
the practical value of A^(fc, i) does not seem compro- 
mised by the minuscule time dependence we have exhib- 
ited. The loop corrections we have discussed can never be 



large (which is the same conclusion reached by Weinberg 
0,Q) because they are suppressed by the quantum grav- 
itational loop counting parameter GH 2 < 1CP 10 . Their 
enhancement by ]n[a(t) / a(tk)] ^ 60 is huge by the stan- 
dards of conventional perturbation theory, and unprece- 
dented in view of its time dependence, but there are sim- 
ply not enough e-foldings of inflation left after first hori- 
zon crossing to overcome the suppression factor for any 
mode whose spatial variation we can now perceive. 

Despite having reached a different conclusion from Sen- 
atore and Zaldarriaga, our results represent no real dis- 
agreement with their analysis. They were uninterested in 
self-interactions from the gravity-innaton system because 
the fixed number of fields in that sector cannot engender 
effects which are enhanced by a potentially large param- 
eter such as Weinberg's J\f. And they dismissed massless 
scalars with nonzero potentials as unnatural. We feel 
it is not reasonable to fine tune the inflaton potential 
V{ip) and then quibble about fine tuning the spectator 
potential U{a). We also thought it worth establishing 
that infrared logarithms do contaminate gauge invariant 
quantum gravity observables such as Aj? [fc , t) because 
the contrary view has been expressed 28j, |29j . 

We close with two thoughts. First, the small infrared 
log corrections to A 2 z (/c,t) might eventually be observ- 
able through 21 centimeter measurements of the matter 
power spectrum out to very large redshifts '30]. This 
would require untangling the primordial signal from late 
time effects, which is very hard but perhaps not impossi- 
ble. It would also require a precise tree order prediction 
from some unique model of inflation. 

Our final comment is that loop corrections to the power 
spectrum are not the best place to study infrared loga- 
rithms because \n[a(t) / a(tk)] cannot exceed about 60 for 
any mode whose spatial variation we now perceive. By 
contrast, there can be spectacular enhancements in quan- 
tities which seem spatially constant, such as the vacuum 
energy [3l| and Newton's constant [j$3, because they 
receive contributions from modes which are still super- 
horizon. For a very long period of inflation perturbation 
theory can even break down, after which reliable com- 
putations would require some nonperturbative resum- 
mation technique. Such a method has been devised by 
Starobinsky [33j, and applied by him and Yokoyama to 
scalar potential models [34| . for which it sums the series 
of leading infrared logarithms [l2[ • Starobinsky's method 
has recently been extended to Yukawa-coupled fermions 
26] and to scalar quantum electrodynamics [35|. It has 
not yet been extended to quantum gravity but there are 
reasons for believing that some version of it can be 18 (31 ] . 
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